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ABSTRACT 

It has been claimed by Taylor et al. that the low-redshift end of the K-z relation for 
radio galaxies is too bright by about half a magnitude due to contributions from the 
obscured quasar nuclei. Such a result has major implications for the use of the K- 
band Hubble diagram in understanding the cosmological evolution of radio galaxies. 
In this paper we present 1-5 /im imaging data of a nearly-complete sample of low- 
redshift radio galaxies; this approach allows us to accurately determine the strengths 
of any unresolved nuclear components in the galaxies. We detect nuclear sources in five 
targets, whose broad-band colours are consistent with reddened quasar spectra. In all 
five cases the ratio of the inferred intrinsic near-infrared luminosity to the narrow-line 
luminosity is typical of quasars. We find a correlation between the inferred nuclear 
extinction and core-to-lobe ratio, which places constraints on the geometry of the 
torus. We find evidence for a shift of the K-z relation to fainter magnitudes, but 
by a much smaller amount (~ 0.1 mag) than Taylor et al. determined. Under the 
assumption that the nuclear sources in radio galaxies have the same intrinsic near- 
infrared spectra as quasars, our multi- wavelength images allow us to limit any possible 
shift to less than 0.3 magnitudes. 
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1 INTRODUCTION 

It has been proposed that powerful FRII (see Fanaroff & 
Riley 1974) radio galaxies and radio loud quasars are intrin- 
sically the same objects, with all observed differences being 
a result of the angle between the line of sight and the radio 
axis (Scheuer 1987; Barthel 1989). While objects classified 
as "narrow-line radio galaxies" do not display obvious broad 
lines in their optical spectra, broad lines have been seen ei- 
ther in the near-infrared (e.g. Hill, Goodrich & DePoy 1996) 
or in polarized light (e.g. Young et al. 1996) . In addition, the 
detection of unresolved nuclei in near-infrared observations 
of radio galaxies (e.g. Djorgovski et al. 1991; Simpson, Ward 
& Wilson 1995), where the optical depth is lower, indicates 
the presence of a quasar-like nucleus in the centres of these 
objects. Since quasars can outshine their host galaxies by 
factors ranging from a few (e.g. Dunlop et al. 1993; Taylor 
et al. 1996, hereafter T96) to more than 100 in extreme cases 
(e.g. PDS 456; Simpson et al. 1999), even if only a few per 
cent of the nuclear light is transmitted it can still make a 
major contribution to the integrated luminosity. 

It has long been known that the brightest radio galaxies, 



i.e. those in the 3CRR catalogue of Laing, Riley & Longair 
(1983), are brighter in the near-infrared by about a magni- 
tude at z « 1 than they are at low redshift. This amount 
of dimming as cosmic time increases is consistent with pas- 
sive evolution of their stellar populations (Lilly & Longair 
1984; Lilly, Longair & Allington-Smith 1985b). While some 
3C radio galaxies at z ~ 1 suffer significant contamination 
from non-stellar radiation (Rawlings et al. 1995; Simpson, 
Rawlings & Lacy 1999), this does not have a significant im- 
pact on the locus of the K—z relation. However, radio galax- 
ies from the fainter 6C and B2 radio surveys are 0.6 mag 
fainter at K than the 3C objects at these redshifts (Eales 
et al. 1997), whereas there is no correlation between ra- 
dio and host galaxy luminosities locally. Furthermore, the 
6C/B2 galaxies lie below the no-evolution curve. 

The possibility that there is significant non-stellar con- 
tamination at low redshift is one that therefore needs to be 
investigated. A study by T96 claimed that the true locus of 
the K—z relation at low redshift may be half a magnitude 
fainter than previously thought due to contamination from 
the non-stellar nucleus. In this paper, we perforin a similar 
analysis to that of T96 on an unbiased sample of radio galax- 
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ies from the catalogue of Laing et al. (1983). In Section 2 we 
describe our sample selection and observational technique. 
In Section 3 we briefly describe the two-dimensional fitting 
procedure used to separate the nuclear and host galaxy prop- 
erties in our J and K images, and present the results of this 
analysis, together with aperture photometry of the longer- 
wavelength images. In Section 4 we use these results to ex- 
amine the nature of the stellar and non-stellar components 
of radio galaxies, and infer the intrinsic luminosities of the 
nuclear sources and the line of sight obscurations towards 
them. In Section 5 we attempt to infer the geometry of this 
obscuring material, and in Section 6 we discuss the impli- 
cations of our analysis on the low-redshift end of the K-z 
relation and compare our results with those of T96. Finally, 
we summarize our results in Section 7. Throughout this pa- 
per we adopt Ho = 50kms -1 Mpc -1 , qo = 0.5, and A = 0. 
Our convention for spectral index, a, is such that the flux 
density at a frequency v, S v oc v~ a . 



2 SAMPLE SELECTION AND OBSERVATIONS 

We selected our targets from the 3CRR catalogue of Laing 
et al. (1983). This catalogue is essentially complete and, by 
virtue of being selected at low frequency, does not preferen- 
tially include sources observed close to the radio axis whose 
fluxes are enhanced by Doppler boosted cores. It is there- 
fore an unbiased sample of radio galaxies, except in the sense 
that it consists of the most radio-luminous objects. We ex- 
cluded objects with known broad Ha, as determined by the 
high-quality spectrophotometric observations of Laing et al. 
(1994, hereafter L94); we did not however exclude 3C 234 
since the broad Ha seen is believed to be entirely scattered 
(Antonucci 1984; Young et al. 1998). This criterion sup- 
plied an upper redshift cutoff of z < 0.43, in order that Ha 
still be present in the optical spectra. From those objects 
visible from the 3.8-m United Kingdom Infrared Telescope 
(UKIRT) in early January, we selected only those galaxies 
with strong optical emission lines (Class A; Hine & Lon- 
gair 1979). The reasons for this were twofold. First, the low- 
ionization ("optically dull", Class B) radio galaxies may not 
be part of the quasar-radio galaxy unification paradigm - 
they are predominantly low-power objects with an FRI ra- 
dio morphology. Secondly, since emission-line luminosity is a 
good indicator of the power of the active nucleus (Rawlings 
& Saunders 1991), it follows that these objects will have less 
luminous nuclei in the near-infrared and are therefore likely 
to be below our detection threshold. We were able to observe 
ten of the eleven objects thus selected, with 3C 123 being 
excluded from our programme solely because it was never 
at a suitable hour angle when the time came to choose the 
next target. We present details of the sample in Table [j], in- 
cluding observed 178 MHz and [O III] A5007 fluxes and rest 
frame luminosities. Although the data of L94 is of a consis- 
tently high quality, their choice of a narrow slit necessarily 
meant a loss of extended emission line flux. This effect is es- 
pecially pronounced in the closest radio galaxies, where the 
line emission can be extended over many arcseconds (see, 
e.g., Baum et al. 1988). We have therefore chosen not to use 
fluxes from L94's data in those cases where a larger aperture 
measurement is available that contains a significantly larger 
flux. We were unable to locate large aperture measurements 



for 3C 42 or 3C 153, although these are the two most distant 
objects in our sample and we can expect the spectroscopic 
slit of L94 to include virtually all the line emission. 

All the images presented here were obtained on UKIRT 
in photometric conditions, using IRCAM3 with a nominal 
scale of 0.286 arcsec pixel -1 . Most observations were made 
on the nights of UT 1996 Jan 3-6, although 3C 33, 3C 153, 
and 3C 171 were observed during an earlier run on UT 1995 
Jan 3-5 as part of a pilot study. Each target was observed 
through the JKL' M filters, with typical on-source exposure 
times of 18 minutes at J and K and 30 minutes at L' and 
M. Most of the radio galaxies are sufficiently compact at J 
and K compared to the size of the array that it was possible 
to use standard jittering techniques and therefore there was 
no need for "off-source" observations to produce a Hatfield 
image. In these cases, each observation consisted of taking 
two sets of nine jittered frames of 60 s duration each (split 
into shorter, background- limited exposures). For those ob- 
jects with z < 0.1, however, nearby regions of sky were used 
to enable flatfielding and sky subtraction. The much higher 
background in the thermal-infrared (L' and M) meant that 
only the central regions of the radio galaxies were detected 
with any significance, and even the closest objects could be 
jittered about the array. However, to permit accurate flat- 
fielding given the rapid background variations, shorter ex- 
posure times were necessary and each observation was split 
into sets of five frames of 30 s (Z/) or 15 s (M), each of which 
was again comprised of a number of shorter, background- 
limited exposures. Due to time constraints, we were unable 
to observe two of our targets (3C 42 and 3C 84) at M. 

Data reduction was performed in the following manner. 
Each observation was split into its sets of nine (at J and K) 
or five (at L' and M) separate exposures, and these were 
dealt with as separate entities. First, an appropriate dark 
frame was subtracted from all the images. Then a flat field 
was produced by median-filtering either the sky frames (if 
the source was at z < 0.1) or the object frames (for cases 
where no off-source observations were made), after scaling 
each frame to have the same median pixel value. Each ob- 
ject frame was then divided by a normalized version of this 
flatfield, and these images were registered, using stars and 
galaxies visible on the individual images where possible, or 
the nominal telescope offsets otherwise. An iterative tech- 
nique was used to determine the DC offsets between regions 
of overlap in separate frames, and the images were averaged. 
At this stage, the separate sets were averaged together to 
produce the final image. The seeing was measured from stars 
in the individual exposures to be in the range 1.0-1.2 arcsec 
FWHM (about 4 pixels) and the registration uncertainties 
did not cause the FWHM of stars in the final coadded images 
to increase by a measurable amount. 

Images of photometric standard stars were taken 
throughout the course of each night. One standard star was 
imaged five times in each filter for each radio galaxy. Mul- 
tiple coadds were used both to improve the signal-to-noise 
ratio, and to ensure that each observation had a similar ex- 
posure time to the individual radio galaxy observations. In 
this way, the point spread function (PSF) of the standard 
star would sample the same longer-term seeing variations 
as the galaxy images. Flux calibration solutions were de- 
termined separately for each night of observation, with the 
r.m.s. dispersion for each night typically being 4% at J and 
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Table 1. The sample of radio galaxies studied. Radio fluxes and spectral indices are from Laing et al. (1983), although the radio fluxes 
have been multiplied by 1.09 to put them on the scale of Baars et al. (1977). The Galactic extinctions used to deredden the [O III] 
luminosities are from NED. References for [O III] flux: YO = Yee & Oke (1978), assuming A5007/A4959 = 3; L = L94 and Laing et al. 
(2000). 
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/5007 


logL 50 07 (W) 


fO TTT1 

[W ill] 


E(B — V) 








(Jy) 




(W Hz- 1 ) 


(lO" 18 Wm" 2 ) 


observed 


dereddened 


ref 




3C 33 


0106+130 


0.0595 


59.3 


0.76 


26.95 


108 


35.20 


35.23 


YO 


0.02 


3C 42 


0125+288 


0.3950 


13.1 


0.73 


27.91 


1.6 


34.96 


35.01 


L 


0.05 


3C 79 


0307+169 


0.2559 


33.2 


0.92 


27.93 


48.5 


36.08 


36.19 


L 


0.10 


3C 84 


0316+413 


0.0172 


66.8 


0.78 


25.93 


544 


34.84 


35.07 


L 


0.17 


3C 98 


0356+102 


0.0306 


51.4 


0.78 


26.32 


76 


34.48 


34.65 


YO 


0.13 


3C 153 


0605+480 


0.2769 


16.7 


0.66 


27.73 


2.2 


34.80 


35.04 


L 


0.23 


3C 171 


0651+542 


0.2384 


21.3 


0.87 


27.68 


21.7 


35.68 


35.74 


L 


0.06 


3C 192 


0802+243 


0.0598 


23.0 


0.79 


26.54 


57 


34.94 


34.98 


YO 


0.03 


3C 223 


0936+361 


0.1368 


16.0 


0.74 


27.10 


43.8 


35.53 


35.53 


L 


0.00 


3C 234 


0958+290 


0.1848 


34.2 


0.86 


27.68 


208 


36.45 


36.47 


YO 


0.01 



Table 2. Photometry in a 3-arcsec aperture. Limits are 3a. 

Galaxy J K 



L' 



M 



3C 33 
3C 42 
3C 79 
3C 84 
3C 98 
3C 153 
3C 171 
3C 192 
3C 223 
3C 234 



14.41 
17.52 
16.37 
13.41 
13.83 
16.71 
16.97 
14.43 
16.02 
15.33 



± 0.04 
± 0.02 
± 0.03 
± 0.03 
± 0.02 
± 0.04 
± 0.04 
± 0.03 
± 0.04 
±0.03 



13.16 
16.06 
14.79 
11.90 
12.60 
15.25 
15.71 
13.46 
14.58 
13.08 



±0.04 
±0.03 
±0.06 
±0.05 
±0.03 
±0.04 
±0.04 
±0.05 
±0.05 
±0.05 



11.89 ± 0.08 

> 15.08 

12.96 ± 0.08 
9.66 ± 0.06 
12.19 ± 0.07 

> 14.97 

> 15.23 
13.39 ± 0.08 
12.91 ± 0.07 
10.50 ± 0.06 



11.17 ± 0.13 
not observed 

> 12.43 

not observed 

> 12.36 

> 12.71 

> 12.44 

> 12.28 
11.79 ± 0.20 

9.91 ± 0.08 



K and 7% at L' and M. Aperture photometry from our im- 
ages agrees with that of Lilly & Longair (1984) and Lilly, 
Longair & Miller (1985a) to within the quoted uncertain- 
ties. In Table ^| we list the photometry measured in a 3- 
arcsec aperture, which is more appropriate for detecting red 
nuclear sources. 



3 SEPARATION OF NUCLEUS AND HOST 
GALAXY 

Thermal emission from the telescope is very strong long- 
ward of 2.5 fim, resulting in a vastly increased background 
and a large drop in sensitivity. In addition, the strength of 
the quasar nucleus relative to the host galaxy increases with 
wavelength due to its red colour. We therefore undertake 
different analyses for the short ( JK) and long (V M) wave- 
length data. 

In our J and K images, the host galaxies are well- 
detected out to large radii and are likely to dominate over 
the nuclei, even in fairly small apertures. Methods which at- 
tempt to estimate the flux of a nuclear source by using aper- 
ture photometry to correct small aperture measurements for 
starlight are prone to overestimate the strength of the source 
(Simpson 1994b) and it is necessary to model the galaxy 
in order to obtain a reliable measurement. At longer wave- 
lengths, the host galaxy is not detected with great signif- 



icance due to the bright thermal background, and simpler 
methods can be employed. 

3.1 Two-dimensional modelling 

To enable the most direct comparison possible of our results 
with those of T96, we have used an almost identical analysis 
technique on our data. We refer the reader to that paper for a 
more detailed description of the method, and for a discussion 
of its advantages over the more common (and simpler) one- 
dimensional analysis (e.g. Simpson et al. 1995). 

The two-dimensional approach involves constructing a 
model image derived from fitting parameters and compar- 
ing this with the actual data, using an error frame to allow 
a quantitative \ 2 minimization. The assumed model was 
an elliptical galaxy obeying a de Vaucouleurs law (de Vau- 
couleurs 1948) and an unresolved nuclear source. There are 
five parameters used to construct the model: 

(i) The luminosity of the nuclear source, L nuc 

(ii) The luminosity of the host galaxy, Lhost 

(iii) The effective (half-light) radius of the host galaxy, r c 

(iv) The position angle of the major axis of the host 
galaxy, 9 

(v) The axial ratio of the host galaxy, a/b 

The pixel location of the centre of the galaxy, (x c ,y c ), 
was determined using iterative centroiding, and construction 
of the model galaxy proceeded as follows. Each pixel in the 
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image was divided into 400 subpixels, and the flux in each 
of these subpixels was computed based on the parameters 
of the galaxy model. This amount of subpixellation ensures 
that the flux in the central pixel is never underestimated 
by more than 0.14%, and any overestimation of the nuclear 
source flux as a result will be negligible. The nuclear source 
was then added by increasing the value of the central pixel 
(x c ,y c ) by the appropriate amount. This nucleus-plus-galaxy 
model was convolved with the adopted PSF to produce a 
simulated image. 

The quality of fit for the simulated image was deter- 
mined by the \ 2 statistic, 

2 ( z{x, y)-C,{x,y) \ 2 

X <r(x,y) J ' 

where z(x,y) and C,(x,y) are the values of the {x,y) pixel 
in the real data and the simulated image respectively, and 
a(x,y) is the uncertainty assigned to this pixel. We con- 
structed an error frame in exactly the same manner as T96, 
although the sampling errors are smaller for our data due 
to the smaller pixel scale. This procedure results in good 
error estimation, and the images of the weighted residuals 
from the best-fit models appeared structureless. Companion 
objects and foreground stars were effectively excluded from 
the fitting procedure at this point by setting the pixel values 
of the error frame to some very large value in regions where 
a contaminating source was visible, thus giving such pixels 
a negligibly small weight in determining the quality of fit. 

In virtually all cases, the radio galaxy occupied only a 
fraction of the total frame. To speed up the computational 
process, we therefore used only a subregion of the image for 
our analysis. We selected a square region which just encom- 
passed the isophote whose intensity was 0.01% of the sky 
brightness, the level at which flatfielding uncertainties are 
expected to become important. Typically such a region was 
about 30 arcsec on a side. The exception was 3C 84, whose 
proximity meant that such an isophote was larger than the 
entire IRCAM3 field of view, and so we used the full array. 
We deviate from the method of T96 in the optimization al- 
gorithm we use. Instead of the Polak-Ribiere minimization 
of T96, we used a simulated annealing algorithm (Press et 
al. 1996). We first tested the fitting procedure on simulated 
galaxies to which noise had been added. After confirming 
that it successfully recovered the input parameters of the 
model, the fitting procedure was performed ten times in J 
and K for each radio galaxy, with two different starting sim- 
plexes for each of the five different PSFs. The results of the 
fit which produced the lowest value of \ 2 were used in later 
analysis. 

We estimate the uncertainties in our fitting parameters 
by considering the variations in the results from the ten sep- 
arate fits. However, not all PSFs were able to produce good 
fits, by which we mean a value of \ 2 close to that of the 
best fit, and we rejected those which did not. As expected, 
only the nuclear luminosity depends strongly on the PSF. 
Typically (for galaxies with a/b > 1.1), the position angle is 
determined to 5° or better, while the uncertainty is as large 
as 10° for 3C 192 (with a/b ~ 1.03, the roundest galaxy in 
our sample). The axial ratio is always determined to better 
than 5%. Although the host galaxy luminosity is generally 
determined to ~ 5% r.m.s., the uncertainty in the effective 
radius is larger, (~ 10%), and the peak surface brightness 



Table 3. Results of the fits to the J and K images. Listed are 
the nuclear and host galaxy (integrated out to infinite radius) 
magnitudes, their ratio in a 12-arcsec aperture, and the effective 
radius (in kpc), axial ratio and position angle of the host galaxy. 



Galaxy m nU c m hoBt f Jnuc (12") r c a/b O xl 

^host 



3C 


33 


./ 


> 18.2 


12.6 


< 


.01 


7.5 


1 


.10 


153 


1 


.78 






A' 


16.6 


11.7 





.02 


5.9 


1 


.09 


158 


2 


.05 


3C 


12 


J 


> 19.0 


16.4 


< 


.13 


15.1 


2. 


.32 


160 


1. 


.58 






K 


> 18.1 


14.9 


< 


.07 


12.9 


1. 


.95 


162 


1. 


.01 


3C 


79 


J 


17.4 


15.5 





.23 


10.4 


1 


.15 


15 


1 


.90 






A 


15.6 


14.2 





.35 


8.0 


1 


.07 


11 


1 


.27 


3C 


84 


J 


15.7 


8.9 


0. 


.02 


24.6 


1. 


.21 


112 


5. 


.53 






A 


12.6 


7.9 


0. 


.15 


25.2 


1 


.26 


114 


8. 


.25 


3C 


98 


J 


> 18.5 


11.9 


< 


.01 


4.2 


1. 


.17 


60 


2 


.72 






A 


> 17.2 


10.8 


< 


.01 


4.5 


1 


.16 


62 


1 


.27 


3C 


153 


J 


> 18.4 


15.4 


< 


.10 


21.5 


1 


.45 


86 





.14 






A 


> 17.5 


14.0 


< 


.07 


21.5 


1. 


.37 


87 


0. 


.22 


3C 


171 


J 


> 18.8 


15.9 


< 


.08 


7.0 


1 


.21 


170 


1 


.24 






A 


> 17.8 


14.5 


< 


.06 


8.6 


1 


.33 


156 


1 


.88 


3C 


192 


J 


> 18.4 


13.0 


< 


.01 


7.8 


1 


.03 


94 


1 


.83 






A 


> 17.4 


12.1 


< 


.01 


7.7 


1. 


.04 


101 


1. 


.65 


3C 


223 


J 


17.5 


14.5 


0. 


.10 


11.1 


1. 


.33 


57 


0. 


.68 






A 


14.9 


13.6 





.50 


13.7 


1 


.42 


59 


1 


.32 


3C 


234 


J 


15.5 


15.4 


1 


.30 


13.7 


1 


.04 


61 


1 


.87 






A 


13.6 


13.8 


1 


.66 


12.2 


1. 


.21 


81 


1. 


.17 



can be uncertain by as much as ~ 20% (with there being 
a strong anti-correlation between effective radius and peak 
surface brightness). The fits we performed in the J and K 
filters for each galaxy were completely independent. How- 
ever, it is expected that some quantities, such as the axial 
ratio and the position angle of the major axis, should be 
very similar in the two bands. The difference between the 
values of these quantities derived for each filter can therefore 
also be taken as a measure of the uncertainty. 

For galaxies with large nuclear-to-host ratios (e.g. 
3C 234) or large angular effective radii (e.g. 3C 84), the 
magnitude of the nuclear source was well-determined. How- 
ever, for other sources it was quite sensitive to the PSF used, 
and in many cases, the best-fit nuclear flux was very faint 
and it became necessary to estimate some limit, fainter than 
which a source would not be considered a true detection. We 
estimate the uncertainty on the nuclear flux (in units of de- 
tector counts) to be a — e^/n, where e is the mean value of 
the error image in a 3-arcsec aperture and n is the number of 
pixels in this aperture. When the flux of the nuclear source 
is less than 5a, we quote an upper limit. These limits vary 
slightly from source to source, depending on such things as 
the seeing, but are typically J « 18.5 and K « 17.5, or a few 
per cent of the integrated host galaxy luminosity. The size 
of the variations in the best fit value of the nuclear flux for 
different PSFs and starting simplexes was broadly consis- 
tent with this level of uncertainty. Contour plots of the data 
and our models can be found in Figs [T|^[lo|- The properties 
of the best-fitting models are listed in Table |j| 

We make a comment regarding the A"-band fit for 3C 98 
listed in Table || This fit has x 2 = 25226/19876, but 
the other nine fits all have x 2 < 26000 yet have nuclear 
sources as bright as K = 14.9. This is in marked contrast to 
the other model galaxies, where large variations in the best- 
fitting nuclear flux compared to the value in Table ^ were 
accompanied by large increases in an d were presumably 
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J (30.75) 3C 33 K (19.75) 




-10 10 -10 10 -10 10 -10 



Figure 1. Images of 3C 33 and the best-fitting model at J (left two panels) and K (right two panels). The number in parenthesis 
indicates the value of the lowest contour plotted (in mag arcsec - 2 ), with successive contours spaced at 0.5 mag intervals. North is up and 
east is to the left. The axes are labelled in arcseconds. These images have not been smoothed. 
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Figure 2. As Fig. m, but for 3C 42. 



J (22.00) 3C 79 K (20.75) 
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Figure 3. As Fig. |, but for 3C 79. 



J (20.25) 3C 84 K (19.25) 




Figure 4. As Fig. |, but for 3C 84. 
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J (20.50) 



3C 98 



K (19.35) 
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Figure 5. As Fig. m, but for 3C 98. 



J (22.00) 



3C 153 
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Figure 6. As Fig. hi but for 3C 153. 



J (22.00) 



3C 171 



K (20.75) 







Figure 7. As Fig. hi but for 3C 171. 



J (21.25) 



3C 192 



K (20.25) 







Figure 8. As Fig. [l| but for 3C 192. 
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Figure 9. As Fig. [l| but for 3C 223. 
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Figure 10. As Fig. ffl, but for 3C 234. 



due to a p oorl y- matched PSF. We shall return to this point 
in Section 4.2. 



Table 4. Inferred magnitudes for the nuclear source in each of 
the radio galaxies. 



3.2 Thermal images 

We were able to detect the radio galaxies at L' in eight 
out of ten cases. We were less successful at M, and only 
detected three sources of the eight which we observed. For 
most of our detections, the flux is dominated by the trans- 
mitted quasar nucleus, although the stellar contribution is 
not negligible. The fact that we measure larger fluxes in 
the larger apertures shows that we have detected starlight 
from the host galaxies, but the signal-to-noise ratio per pixel 
is far too low to allow the two-dimensional modelling used 
for the shorter wavelength images. Indeed, for most galaxies 
the extended emission is not bright enough to even permit 
a meaningful measurement of the flux in a circumnuclear 
annulus. We must therefore estimate the stellar flux by ex- 
trapolating from our shorter wavelength images. 

The stellar populations of early-type galaxies are domi- 
nated in the near-infrared by stars on the red giant branch. 
We therefore use the observed spectrum of the M2II-III star 
(3 Peg (Strecker, Erickson & Witteborn 1979) as a model for 
the stellar populations of our galaxies. The colours of (3 Peg 
are known to match those of early-type galaxies, and the 
star's spectrum has been used to compute if-corrections for 
elliptical galaxies (see Longmore & Sharpies 1982). We cor- 
rect the observed 3-arcsec K magnitudes for the presence 
of the nuclear source listed in Table ^ to determine the flux 
from stars only in this aperture, and use the (3 Peg spectrum 
to estimate the stellar contribution to the flux at L' and M. 



Galaxy 




./ 




A' 




V 
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3C 33 
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18.20 




16.57 




12.32 




11.33 


3C 42 
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19.00 


> 


18.10 


> 


15.82 






3C 79 




17.37 




15.89 




13.11 
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12.54 


3C 84 




15.70 




12.59 




9.73 






3C 98 
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18.50 


> 


17.20 
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13.10 
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13.46 


3C 153 
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18.40 


> 


17.50 


> 


16.58 


> 


12.86 


3C 171 


> 


18.80 


> 


17.80 


> 


16.58 


> 


12.51 


3C 192 


> 


18.40 


> 


17.40 


> 


14.99 


> 


12.63 


3C 223 




17.52 




14.91 




12.98 




11.81 


3C 234 




15.54 




13.65 




10.56 




9.94 



If there is an excess at greater than 3<r significance in the 
3-arcsec aperture, we attribute it to the quasar nucleus. The 
results of this analysis are presented in Table ^| 

Although 3C 98 and 3C 192 are both detected at L', 
the emission we observe is purely stellar in origin, as both 
sources have 3-arcsec aperture K — L' colours consistent 
with the (3 Peg spectrum, although the colour of 3C 98 is 
redder than f3 Peg's and only marginally consistent with it. 
However, Lilly et al. (1985a) showed that even in a much 
larger aperture (7.5arcsec), 3C 98 is redder than a typical 
early-type galaxy, which rules out a compact red source as an 
explanation for the colour. We note that the L' — M < —0.17 
limit is also consistent with starlight, but is far too blue for 
a quasar (corresponding to a power law index a < —2.5). 
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Figure 11. Two colour diagram for the ten radio galaxies. 
Colours have been measured in a 3-arcsec aperture. 



3.3 Reality of our detections 

If we examine Fig. [nj which presents the 3-arcsec aper- 
ture colours of the ten radio galaxies from Table |^, the 
galaxies split naturally into two groups: those with K — L' 
colours consistent with late-type stellar populations, and 
those which are much redder. There is also a weak, yet highly 
significant, correlation between the two colours, as expected 
from a red nuclear component which becomes increasi ngly 



3.2 



significant at K. The simple colour analysis of Section 
will obviously infer a nuclear source for galaxies in the sec- 
ond group only, but it is important to note that these are also 
the only galaxies for which nuclear sources were detected in 
Section to. All nuclear sources which are reliably detected 
at K in our 2D analysis should be detectable in our L' im- 
ages due to their red colour; this is indeed true. Similarly, 
any real nuclear sources we detect at J should also be seen 
at K (and therefore at L' too), and this is also the case. 
A similar argument does not apply to our M-band data, 
since the reduction in sensitivity over L' (~ 3 mag) is much 
greater than can be compensated for by the red colour of 
the quasar nucleus (Ll — M = 0.76 for a power law a = 1.3; 
Neugebauer et al. 1987). We are therefore convinced that 
we have not made any spurious detection, and that we have 
not failed to detect any sources which are in reality brighter 
than our quoted limits. 



4 LUMINOSITY AND EXTINCTION 
ESTIMATES 

4.1 Galaxies with detected nuclei 

If the standard interpretation of our nuclear sources as dust- 
obscured quasars is correct, they should have the colours of 
a quasar, modified by some amount of foreground reddening. 
We therefore perform linear regression, accounting for upper 
limits in the manner of Isobe, Feigelson & Nelson (1986), on 
the data of Section ^j. This analysis produces values for the 
intrinsic luminosity and extinction of each nuclear source. 
We adopt our own parametrization of the near-infrared in- 
terstellar extinction law, obtained by fitting a second-order 
polynomial to the data of Rieke & Lebofsky (1985). This 
parametrization, 



Table 5. Foreground extinction and unobscured rest-frame 1 /im 
luminosity for the obscured quasar nuclei detected in our analy- 
sis. The dereddened [O III] A5007 emission line luminosity from 
Table |l| is also listed. 



Galaxy 


Ay (mag) 


logii Mm (WHz- 1 ) 1 


Ogi5007 (W) 


3C 33 


29.5 ±3.5 


22.49 ± 0.11 


35.23 


3C 79 


1.0 ±2.5 


23.07 ±0.25 


36.19 


3C 84 


9.4± 1.6 


22.02 ± 0.12 


35.07 


3C 223 


3.9 ±0.6 


22.72 ± 0.05 


35.53 


3C 234 


4.0 ±1.9 


23.83 ± 0.15 


36.47 



0.243A" 2 + 0.184A" 1 - 0.022, 



by virtue of ignoring the optical data, provides a rather bet- 
ter fit to the longest wavelengths than do those of Howarth 
(1983) and Cardelli, Clayton & Mathis (1989). The intrinsic 
quasar spectrum is assumed to be that of an a = 1.3 power 
law (Neugebauer et al. 1987). Table |E] lists the extinctions 
and intrinsic luminosities we derive for the obscured quasars 
in the five radio galaxies with detected nuclear sources. 
These data are presented graphically in Fig. For 3C 223 
and 3C 234, where the nucleus was detected at all four wave- 
lengths, we also performed regression with the spectral index 
as an additional free parameter. In both cases the spectral 
index was found to be close to 1.3 and the results were not 
changed significantly. 



4.2 Galaxies without detected nuclei 

The five galaxies with detected nuclei will obviously prefer- 
entially include those objects with low nuclear extinction. If 
we wish to reach some conclusions about the amount of ob- 
scuration in a 'typical' radio galaxy, we must make estimates 
of the extinction in the remaining five members of the sam- 
ple. This requires us to calculate the intrinsic near-infrared 
luminosity of the quasars in these objects, which we do in 
the following manner. 

We first estimate the unobscured optical luminosity of 
the quasar using the [O III] A5007 emission line. Simpson 
(1998a) has argued that this line is a good indicator of the 
intrinsic luminosity of radio galaxies and lobe-dominated 
quasars (such as those selected at low radio frequency). It 
has also been shown that the [O III] luminosities of lobe- 
and core-dominated quasars are the same for a given ex- 
tended radio luminosity (Jackson et al. 1989; Jackson & 
Browne 1991; Corbin 1997), indicating that the line is emit- 
ted isotropically. 

There is no Baldwin effect in the equivalent width of 
[O III] in lobe-dominated quasars (Jackson & Browne 1991). 
Only when core-dominated objects are considered does such 
an effect appear. This can be understood as the result of an 
increase in the continuum luminosity at small polar angles 
due to an anisotropic component. However, since the near- 
infrared emission is believed to be produced by dust heated 
by the optical-ultraviolet radiation and located close to the 
equatorial plane of the quasar, it is the luminosity at large 
viewing angles, such as is seen in lobe-dominated quasars, 
which is of interest to us. 

By virtue of the above arguments, we can use the 
optical-infrared properties of lobe-dominated quasars as 
benchmarks for the properties of the quasars within our ra- 
dio galaxies. The mean [O III] equivalent width for the R < 1 
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Figure 12. Rest-frame nuclear luminosities and best-fit reddened 
quasar models for the five radio galaxies in which a nuclear source 
was detected. For clarity some of the galaxies have been shifted 
upwards in the y-direction by the number of dex given in paren- 
theses. 

quasars from Jackson & Browne (1991) is logEWsoo7 (A) = 
1.66 ± 0.32, or 46 A with a factor of 2 scatter. This is larger 
that the equivalent width of 24 A determined by Miller et 
al. (1992) for low redshift (z < 0.5) quasars selected from 
the Palomar-Green Bright Quasar Sample (BQS; Schmidt 
& Green 1983), because optically selected samples like the 
BQS preferentially contain objects oriented close to the line 
of sight where the optical continuum is boosted. This is sup- 
ported by the large values of the 5-GHz core-to-lobe ratio, 
R, displayed by the radio-loud quasars in the low redshift 
BQS (Miller, Rawlings & Saunders 1993). 

Adopting 46 A as the rest-frame equivalent width (with 
respect to the non-stellar continuum) of the [O III] A5007 
emission line, and assuming a = 0.4 ± 0.4 as representative 
of the continuum between 5007 A and 1 /im (Neugebauer et 
al. 1987) , we find that the last two columns of Table |H| should 
be related by 

log Li Mm - log Lsoot = 12.62 ± 0.34 , 

which is indeed the case, and indicates the validity of our 
method. 
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Figure 13. Inferred nuclear extinction for the galaxies in our 
sample, plotted against the core-to-lobe ratio, R. An anti- 
correlation is readily apparent. 



Applying this technique to the five radio galaxies with- 
out detected nuclear sources, we derive lower limits to the 
nuclear extinction. The non-detections of compact sources 
in 3C 42 and 3C 153 do not rule out unobscured nuclei with 
any significance, since both these sources have faint emission 
line fluxes and therefore are expected to have quite faint nu- 
clei. For the other three galaxies we determine the following 
3cr limits: A v > 74 m for 3C 98; A v > 37 m for 3C 171; and 
A v > 45 m for 3C 192. Note that these limits are derived 
from the non-detections at M-band (where our sensitivity 
to highly obscured nuclei is greatest) and are therefore in- 
dependent of the J and 7\-band fits. The high extinction 
towards the nucleus of 3C 98 therefore clearly precludes a 
JT-band nucleus brighter than the limit we quote in Table [| 



5 THE GEOMETRY OF THE OBSCURING 
MATERIAL 

The core dominance parameter, R (also called the core-to- 
lobe ratio), is a widely-used orientation indicator (Orr & 
Browne 1982) and we are therefore capable of determining 
the nuclear extinction to our radio galaxies as a function of 
viewing angle. Simpson (1994a, 1996) proposed analysing a 
large sample of radio galaxies in this way to distinguish be- 
tween different geometries for the obscuring material, which 
partly provided the motivation for the study described here. 
The value of R is usually determined at 5 GHz, and the prop- 
erties of our sources at this frequency are listed in Table Al 
of the Appendix. We omit 3C 84 since it does not possess 
the characteristic FR II radio morphology. The values of the 
nuclear extinction derived earlier are plotted against log R 
in Fig. || 

Using a generalized Kendall's rank correlation statistic 
(Isobe et al. 1986), we can reject the hypothesis that Av is 
uncorrelated with R at better than 95% confidence, even if 
3C 234 is excluded from the analysis. The nuclear extinction 
therefore increases with viewing angle (i.e. as R decreases) 
and suggests a flattened distribution for the obscuring ma- 
terial. 

Ignoring density inhomogeneities, the obscuring mate- 
rial's geometry is defined by Av(0), which we can derive 
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Table 6. Core dominance parameters and viewing angles with 
respect to the radio axis for the radio galaxies in our sample. 



Galaxy log R 



(°) 



3C 33 
3C 42 
3C 79 
3C 98 
3C 153 
3C 171 
3C 192 
3C 223 
3C 234 



-2.33 
-2.69 
-2.19 
-2.75 
-3.52 
-2.82 
-2.55 
-2.20 
-1.29 



57; 
64; 
53; 
65; 

> 54 
67; 

62; 
54; 

33 



-16 
-12 
-11 

-13 

-18 

-11 

-12 

-14 

(2o) 

-10 

-16 

-13 

-14 

-17 

-12 

-Ifi 

-7 



from Av(R) if we can transform the observed core domi- 
nance parameter, R, to a viewing angle, 8. These two quan- 
tities are related by 



R = -Rt{1 - /3cos6>)~ 



+ 2-Rt(1 + /3cos( 



(Orr & Browne 1982), where R T = R(90°) and /3 is the 
velocity of the beamed core component relative to the speed 
of light. We adopt (3 = 0.99 (Orr & Browne 1982), although 
its exact value is unimportant except at small values of 9, 
which are not relevant to our study of radio galaxies. 

To convert from 71 to 8, it is necessary to know the value 
of Rt- Orr & Browne determined Rt = 0.024 from a sample 
of 3CR quasars, but they assumed that these objects were 
randomly oriented with respect to the line of sight. Since it is 
now believed that quasars are preferentially pointed towards 
us, Orr & Browne's value of Rt will be an overestimate. We 
therefore redetermine Rt using both quasars and FRII ra- 
dio galaxies from Laing et al. (1983) with z < 0.43, although 
we exclude 3C 236 since it possesses a steep-spectrum radio 
core. As described in the Appendix, a good fit to the data 
is obtained with P(logi? T ) = N(-2.54, 0.51). This is con- 
sistent with the distribution derived by Simpson from the 
radio galaxies alone, and with estimates made from other 
samples of radio galaxies (e.g. Morganti et al. 1997). 

Since we have only determined the a priori probability 
distribution for Rt, and not its actual value for each ra- 
dio galaxy, the observed core dominance parameter, R, is 
transformed into a probability distribution for the viewing 
angle, 8. This transformation is described in the Appendix. 
We list the viewing angles for the radio galaxies in Table ^ 
and plot the cross-section through the torus in Fig. ju| It 
is immediately clear from this analysis that our inability 
to constrain the angle at which we are viewing each object 
is a major limitation to the effectiveness of this technique. 
Doppler boosting only produces a fourfold increase in the 
observed luminosity of the central component for 8 « 50° 
(the viewing angle which is believed to separate quasars from 
radio galaxies at low redshift; L94, Simpson 1998a, and see 
also Table |^) , barely enough to produce a signal over the fac- 
tor of three scatter in Rt- Since core dominance is believed 
to be the best available orientation indicator, it is clear that 
a much larger sample of radio galaxies will be needed to 
produce a useful torus map. 
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Figure 14. A cross section through the torus perpendicular to 
the equatorial plane. The axes are in units of magnitudes of visual 
extinction. The dotted lines describe the polar coordinate frame 
and are there to guide the eye. 



6 THE K-Z RELATION AT LOW REDSHIFT 

T96 claimed that the K magnitudes of their radio galaxy 
hosts are 0.5 magnitudes fainter than the currently- assumed 
K—z relation, with the discrepancy being caused by the con- 
tribution from the quasar nucleus. They detected nuclear 
sources in all twelve of their galaxies, with a mean nuclear- 
to-host ratio (Lnuc/Lhost) = 0.80. Our study, on the other 
hand, failed to detect a quasar contribution in five of our 
objects, despite a similarly faint detection threshold. What 
effect does our result have on the K—z locus? We find the 
best least-squares fit to our data that can be obtained by 
applying an offset, A, to the K—z relation of Lilly et al. 
(1985b), 



K = 49.87 - a/1092.7 - 259.74 logz + A. 

For the raw data, the best offset is effectively zero (A = 
—0.01), as expected given the accuracy of our photometry 
and the unbiased nature of our sample. If we consider instead 
just the host magnitudes of Table H, we find A = 0.10, sup- 
porting a shift of the K-z relation to fainter magnitudes, but 
by a much smaller amount than that proposed by T96. In 
addition, this should probably be considered an upper limit 
to any shift required in the K—z relation, since sources from 
fainter radio surveys (which lie on same K—z relation due 
to the lack of a correlation between radio and host galaxy 
luminosity and low redshift) will have intrinsically fainter 
infrared nuclei, which should therefore contribute less flux. 
We display our results in Fig. [l5| together with the host 
galaxy magnitudes derived by T96, which are clearly fainter 
on average. We investigate this discrepancy, first by address- 
ing the difference in the detection rates for nuclear sources 
between our study and that of T96, and then by looking at 
the difference in host galaxy magnitudes. We assume that 
the nuclear sources we and T96 have detected are in fact 
active nuclei, based on their luminosities and broad-band 
SEDs (T96 make a similar claim) . If instead they are merely 
stellar cusps (albeit ones with unusual colours), then there 
is obviously no reason to shift the K—z relation at all, since 
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Figure 15. The if -band Hubble diagram for low redshift radio 
galaxies. Open stars are the observed magnitudes for our 3CRR 
radio galaxies, filled stars are the model host galaxy magnitudes 
for these objects, and crosses and open squares are the observed 
and model host galaxy magnitudes for T96's radio galaxies (all 
measured or synthesized in 12-arcsec apertures). The dotted line 
is the K—z relation from Lilly et al. (1985b), and the solid line is 
this relation shifted by 0.10 mag to provide the best least-squares 
fit to our model magnitudes. 

the actual form of the radial surface brightness profile is 
unimportant. 

6.1 Detection efficiency 

While our sample is effectively complete (only the ran- 
dom exclusion of 3C 123 prevents it from being truly so), 
the twelve objects studied by T96 were selected to match 
previously-constructed samples of radio-loud and radio- 
quiet quasars. One criterion was a match in radio spectral in- 
dex, and T96's sample therefore included a number of radio 
galaxies with spectra much flatter than the typical a « 0.7. 
There are five such fiat-spectrum radio galaxies in the sam- 
ple of T96, and they are therefore strongly overrepresented 
compared to the general low frequency selected radio galaxy 
population (there are none in our sample). Flat integrated 
radio spectra are produced when the flat spectrum radio core 
dominates over the steep spectrum lobes in the frequency 
region of interest, and consequently such sources have high 
values of the core-to-lobe parameter, R. Fig. |l3| showed the 
strong anti-correlation between Av and R, which we inter- 
preted in terms of the geometry of the obscuring material. 
We should therefore expect the nuclei of T96's radio galaxies 
to be less heavily obscured than those of our sample mem- 
bers. This should result in brighter (and hence more readily 
detectable) nuclei and a tendency for their sources to lie be- 
low (i.e. be brighter than) the K-z relation. Indeed, an offset 
of 0.14 mag to brighter magnitudes needs to be applied to 
the relation of Lilly et al. (1995b) in order to obtain the best 
fit to their raw data. 

The single-filter imaging of T96 does not permit the de- 
tailed method of extinction determination we employed in 
Section ^l]. Instead we estimate the average extinction for 
T96's radio galaxy sample as a whole, by comparing the lu- 
minosities of the nuclei detected in the radio galaxies with 



those detected in the radio-loud quasars. Due to the way the 
samples were matched, the intrinsic luminosities of the ac- 
tive nuclei should be the same. Their Tables 4 and 5 indicate 
that the mean observed luminosity of the radio galaxy nuclei 
is 1.77 mag fainter than that of the radio- loud quasar nu- 
clei. A similar result is found when the Lnuc/^host ratios are 
compared. For a typical redshift z « 0.2, this corresponds 
to Av ~ 13 m . By contrast, our sample has an average ex- 
tinction A v > 23 m . 

Our estimate appears to be more in line with the ad- 
mittedly very limited previous study in this field. Simpson 
(1994a) studied a small, but complete, sample of nearby 
radio galaxies from the MS4 survey (Burgess & Hunstead 
1994) which included three Class A (Hine & Longair 1979) 
objects. A nuclear source was detected in one object, which 
further study revealed to suffer an extinction Ay = 30 ± 4 m 
(PKS 0634-205; Simpson et al. 1995). The lower limits de- 
termined for the other two Class A sources implied an av- 
erage Av > 21 m . The powerful radio galaxy Cygnus A has 
also had its quasar nucleus detected in the near-infrared 
(Djorgovski et al. 1991), and the most recent extinction es- 
timates indicate very heavy obscuration (Av ~ 150 m ; Simp- 
son 1994a, Ward 1996). 

We also note that two of T96's radio galaxies, 3C 219 
(L94; Hill et al. 1996) and 3C 287.1 (Eracleous & Halpern 
1994) clearly show broad Ha in their optical spectra. It 
is therefore to be expected that these two objects possess 
bright nuclear sources in the near-infrared. We conclude that 
T96's sample is predisposed towards sources with relatively 
low nuclear obscuration, and believe that our sample, being 
effectively complete and free from selection biases, is more 
representative of the radio galaxy population at large. 

6.2 Host galaxy magnitudes 

Even though T96's sample of radio galaxies appears to be 
biased towards objects with low nuclear obscuration, the 
host galaxies should not be affected. Our multi-wavelength 
images allow us to compute robust upper limits to the non- 
stellar flux at K (and hence an upper limit to the shift of 
the K-z relation required) in a relatively simple manner 
which does not involve a "black box" routine such as the 
two-dimensional fitting procedure. First, we assume that all 
the L' flux in a 3-arcsec aperture arises from the nucleus - 
since there must be some starlight present, this will over- 
estimate the true nuclear flux. In the three instances where 
no detection was made at L' , we use the 3cr upper limit 
(again, overestimating the true flux). We then determine an 
upper limit to the K-band flux by assuming that the spec- 
trum is an unreddened a = 1.3 power law. Although there 
will be variations in the near-infrared spectral indices of the 
objects in our sample, these should average out among the 
ten galaxies. We compute the faintest possible host galaxy 
K-band magnitude by subtracting this nuclear flux from the 
observed 12-arcsec aperture flux. When this analysis is ap- 
plied to 3C 234, the maximum nuclear strength we derive 
exceeds the observed 12-arcsec magnitude, which indicates 
the conservative nature of our approach. Obviously, we can- 
not merely exclude 3C 234 from the analysis, since it is the 
most strongly contaminated source in our sample, so we as- 
sume the host has K = 14.2, as determined by our model- 
fitting (T96 obtain a very similar result). Even with these 
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gross overestimates of the nuclear fluxes (and hence under- 
estimates of the host galaxy brightnesses), the offset from 
the K-z relation is only A = 0.35, somewhat less than the 
0.5 mag found by T96. More realistic estimates, such as as- 
suming that the nucleus suffers at least 5 mag of extinction 
(just enough to obscure broad Ha), imply A < 0.27. The K— 
z relation followed by the host galaxies in our sample cannot 
therefore be the same as that found by T96 for the galaxies 
in their sample, unless the nuclei of our radio galaxies have 
intrinsic near-infrared spectral indices of a < 0.2. This is 
obviously incompatible with unification models which claim 
that the nuclei of radio galaxies are identical to those of 
quasars. 

How then did T96 find such a large shift in the K-z 
relation? We first rule out the possibility that the galaxies 
in T96's sample are genuinely fainter than ours. Although 
T96's sample is generally of lower radio luminosity than 
ours, Hill & Lilly (1991) have shown that there is no differ- 
ence in the optical magnitudes of radio galaxies over three 
orders of magnitude in radio luminosity, and the optical- 
infrared colours of the two samples cannot differ by half a 
magnitude. In addition, Eales et al. (1997) have shown that 
the K magnitudes of radio galaxies from the 3C and 6C 
samples (6C is selected at a ~ 5 times fainter flux level) 
are similar at low redshift. Finally, if we exclude those ra- 
dio galaxies from T96's sample which do not meet the 3CRR 
flux limit, their result is virtually unchanged. Similarly, if we 
consider only FRII galaxies from T96, or only those galaxies 
with a steep radio spectrum, we find no significant change 
to their result. 

We consider two likely explanations for the different K— 
z relations inferred. We first consider what effect the small 
field of view (~ 35") of IRCAM1 might have had on T96's 
modelling. While the much larger field of view of IRCAM3 
provides a significant amount of blank sky in our images 
and therefore allows an accurate determination of the back- 
ground level, this was not the case with IRCAM1. As T96 
explain, they may have overestimated the sky level in their 
images, and we discuss here why this can lead to an under- 
estimation of the host galaxy brightness. T96's radio galax- 
ies have a typical effective radius of a few arcseconds. At 
the edge of their field of view, a de Vaucouleurs profile does 
not have a negligible surface brightness, but rather will have 
/i « 23 mag arcsec -2 (for a typical ~ 20.5 and r c w 6"). 
We have investigated the effects of this on modified K-btxnd 
images of some of our galaxies, by subtracting a constant 
value from all the pixels to simulate overestimation of the 
sky level. 

We find that for small deviations from the true value, 
the effective radius of the best-fitting model does not change 
significantly, but the flux normalization of the stellar com- 
ponent is decreased. This normalization is primarily deter- 
mined at r ~ r c , since the signal-to-noise ratio is much lower 
at larger radii, while the sampling errors can be large near 
the centre, and the flux of the nuclear source can be varied 
to compensate for deficiencies in the quality of the fit. For a 
typical galaxy with an error in background determination of 
the level indicated above, we find an increase of ~ 20-40% 
increase in the flux attributed to the nucleus. The value of 
the \ 2 statistic also increases, and it is therefore worth con- 
sidering whether more accurate results might be obtained if 



Table 7. Direct comparison between our results and those of the 
elliptical galaxy fits of Taylor et al. (1996). Host galaxy magni- 
tudes are within a simulated 12-arcsec aperture. 



Object 




^nuc 


-K"host 


Lnuc 
-khost 


r e 


a/b 


e 


3C 79 


Us 


15.9 


14.4 


0.35 


8.0 


1.07 


ii 




T96 


15.2 


14.6 


0.57 


16.6 


1.16 


7 


3C 234 


Us 


13.6 


14.2 


1.66 


12.2 


1.21 


81 




T96 


14.0 


14.3 


1.41 


13.0 


1.22 


55 



the sky level is included as an extra fitting parameter in the 
modelling process. 

There is also evidence that T96 may have systemati- 
cally overestimated the scale lengths of their galaxies. Opti- 
cal imaging with the Hubble Space Telescope (Dunlop et al. 
2000) produced a median scale length of 12 kpc, compared to 
20kpc determined by T96. The effect of overestimating the 
r e is also to lower the contribution from the host galaxy in 
the central regions, and therefore to overestimate the con- 
tribution from a nuclear point source. A systematic error 
of ~ +70% in the scale lengths determined by T96 would 
serve to produce central surface brightnesses ~ 1 mag too 
faint (again, the fitting procedure will tend to match the 
model to the data at r w r c ) and the model flux within a 
12" aperture will be underestimated by ~ 0.4-0.7 mag, de- 
pending on the size of the galaxy. Since this is very much in 
line with the shift determined by T96, it would seem to be 
a likely explanation for their result. 

We can make a direct comparison between our results 
and those of T96 for the two objects which are common 
to both samples. We present the results of our 7f-band fits 
and those of T96 in Table @. The agreement for 3C 234 
is excellent, but it is less good for 3C 79. Optical imag- 
ing with the Hubble Space Telescope (Dunlop et al. 2000) 
produces an effective radius r c — 9.4 kpc, somewhere be- 
tween the values we determine from our J and Jf-band im- 
ages, and clearly indicating that T96 overestimated the true 
value. A visual inspection of their Fig. A2 appears to indi- 
cate that the larger effective radius is influenced by the com- 
panion source approximately 6" north of the radio galaxy. 
Alternatively, there may have been additional scattered light 
within IRCAM1 which would artificially enhance the surface 
brightness at large radial distances and hence cause the scale 
length to be incorrectly overdetermined, although the good 
agreement for 3C 234 (where the bright nuclear point source 
would be expected to produce a large effect) seems to argue 
against this explanation. 

It is quite possible that different effects are at work for 
different objects. Since we can think of no likely effect which 
would serve to overestimate the host galaxy brightnesses, 
we can expect a significant systematic error in the results of 
T96, even if we cannot be certain what the dominant cause 
of this error might be. 



7 SUMMARY 

We have analysed near-infrared JKL'M images of an ef- 
fectively complete sample of ten 3CRR radio galaxies with 
strong emission lines, and have separated them into their nu- 
clear and stellar components. We find that the colours of the 
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nuclear components are well-modelled by a reddened power- 
law, and the derived unobscured luminosities are in excellent 
agreement with the values predicted from the emission-line 
luminosities. 

We find evidence that the nuclear extinction increases 
with viewing angle, suggesting a flattened distribution for 
the obscuring material. The contributions from the nuclear 
components make the model host galaxy magnitudes fainter 
than the observed K magnitudes, and support a shift of 
the K—z relation to fainter magnitudes by 0.1 mag for low- 
redshift 3CRR galaxies. Using our observed L' magnitudes 
and the assumption that the nuclear sources in radio galaxies 
are intrinsically no bluer than quasar nuclei, we can rule 
out a shift of more than 0.3 mag, independent of our two- 
dimensional model-fitting procedure. 

We have investigated possible causes of the different re- 
sults obtained by Taylor et al. (1996) and ourselves. While 
we believe that T96's sample is biased towards sources with 
low nuclear obscuration and therefore a relatively large nu- 
clear fraction, this should not affect their host galaxy mag- 
nitudes. We believe that their fainter magnitudes are caused 
by overestimating the sky level in their images and/or the 
effective radii of their galaxies. We conclude that the mag- 
nitudes of low-redshift radio galaxy hosts have been system- 
atically overestimated by no more than about 0.1 mag as a 
result of non-stellar radiation from the obscured quasar nu- 
clei, and that conclusions drawn so far concerning the cosmic 
evolution of radio galaxies (e.g. Lilly & Longair 1984; Eales 
& Rawlings 1996) are therefore reliable. 
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APPENDIX A: CORE DOMINANCE AND 
VIEWING ANGLE 

This Appendix describes how we have used the observed core 
dominance parameter of a radio galaxy to infer the viewing 
angle between its radio axis an d o ur line of sight. To do 
this, we have compiled in Table Al radio data for all FRII 
sources with z < 0.43 in Laing et al. (1983). In the interests 
of keeping a homogeneous dataset, we use the 5- GHz flux 
measurements and spectral indices between 2.7 and 5 GHz 
from Pauliny-Toth & Kellerman (1968). Although these are 
not necessarily the most accurate data available, this is not 
a concern as our assumption of a common spectral index 
a — for the core components (which may be somewhat 
variable) will introduce some errors in our determination of 
the rest-frame 5-GHz core dominance parameter. For the 
four non-3C sources which were not observed by Pauliny- 
Toth & Kellerman, we have opted to use the data from 
Herbig & Readhead (1992). These authors perform fits to 
the radio spectra of sources, which should also help to ho- 
mogenize the data they compile from the literature. Her- 
big & Readhead list luminosities and spectral indices at a 
rest-frame frequency of 2.5 GHz, which we use to compute 
a rest-frame 5-GHz flux. However, for consistency in Ta- 
we list synthetic observed-frame 5-GHz fluxes deter- 



Hine & Longair (1979) Class B radio galaxies in this sam- 
ple to improve the statistics (the published spectroscopy of 
the late-RA sources is too poor to make completely accurate 
classifications and we would therefore have to exclude all of 
them), their R distribution is not significantly different from 
that of the Class A objects (a result also found by Morganti 
et al. 1997). 

The first step in our procedure is to determine the dis- 
tribution of the transverse core dominance parameter, Rt- 
We construct synthetic samples of 5000 radio sources which 
are randomly oriented with respect to the line of sight, and 
have the logarithms of their transverse core dominance pa- 
rameters, logi?T, drawn from a Normal distribution, whose 
mean and standard deviation we vary to determine the best 
fit to the observed data. We quantify the goodness-of-fit us- 
ing the Peto-Prentice generalized Wilcoxon test (Latta 1981; 
Feigelson & Nelson 1985), as employed in the ASURV pack- 
age within iraf/stsdas. Of the various tests Feigelson & 
Nelson list to compare two samples of censored data, this 
is believed to be the least affected by differences in sample 
size and censoring patterns, which is important here since 
our synthetic sample is much larger than the real dataset 
and suffers no censoring. The best fit we find is when the 
transverse core dominance parameter is drawn from the dis- 
tribution N(— 2.54, 0.51), with the probability that the two 
distributions are drawn from the same population being 
99.9%. Both versions of the Gehan test (also available in 
ASURv) produce the same result, so it is not sensitive to 
the specific test used. Note that the underlying distribution, 
log7? T G N(-2.43, 0.50), reported by Simpson (1998b) was 
computed using a few older core flux measurements (gener- 
ally upper limits which have been superseded by more sensi- 
tive limits or faint detections) , and with a different goodness- 
of-fit statistic. However, there is only a slight difference be- 
tween the two distributions, which does not affect the results 
in that paper. 

Having derived a distribution function for logi?T, it be- 
comes possible to infer the probability distribution of view- 
ing angles, P(6\R). Bayes' Theorem informs us that 

P{6\R) oc P{R\6)P{6). 

Since R and 6 are directly related by Rt, we can write the 
first term on the right hand side as an a priori probability 
that Rt lies in the appropriate range to produce the ob- 
served value of R given that the viewing angle lies in the 
interval between 9 and 9 + 59. If we define x = f3 cos 9, it is 



trivial to show that R/Rt = (1 + x )(1 



Hence, 



P(7?|6?)=P(lo gj RT = lo gj R 1+ , 

and, straightforwardly, 
d(logi? T ) I3sin9 2x(3 + x 2 ) 



T - x 2 ) 2 dlog Rt 



dO 



A9 



In 10 



1 



ble AT 



mined from their fits. Note that although we have included 



Finally, we need to consider the a priori distribution 
of viewing angles, P{9). A natural choice, motivated by the 
radio galaxy-quasar unification scenarios, would be to as- 
sume that P(9) oc sm9 for angles greater than some critical 
angle, 9 C , and zero for angles less than this. In practice, the 
actual choice of 9 C has little effect on the results, except 
for sources with bright radio cores, such as 3C 234, when 
the probability distribution of angles piles up at the critical 
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Table Al. Properties of 3CRR FRII radio sources with z < 
0.43. The total 5-GHz fluxes and spectral indices between 2.7 and 
5 GHz are from Pauliny-Toth & Kellerman (1968) except where 
an asterisk (*) denotes synthesized data from the spectral fits 
of Herbig & Readhead (1992). We preferentially use 5-GHz core 
flux measurements, where available, and assume a flat spectrum 
to convert measurements at other frequencies. The last column 
identifies which paper in the reference list was used to determine 
the core flux. Fluxes cited as from LP (Leahy & Perley 1991) are 
estimated from their radio maps. 



Name 



otot 
°5GHz 

(Jy) 



gcorc 

(mjy) 



logR 



Core 
rcf 



3C 


16 


0.4050 


0.51 


0.99 


< 0.5 




—3.15 


LP 


3C 


20 


0.1740 


4.18 


0.70 


2.6 




—3.25 


FBP 


3C 


33 




5.03 


0.57 


24.0 




—2.33 


HM 


3C 


33.1 


0.1810 


0.86 


1.13 


15.0 




— 1.83 


vBJ 


3C 


35 


0.0670 


0.59 


1.31 


14.0 




— 1.65 


vBJ 


3C 


12 


3950 


0.84 


1.03 


2.4 




—2.69 


FBP 


3C 


47 


0.4250 


1.10 


1.01 


80.0 




— 1.27 


PH 


3C 


61.1 


0.1860 


1.91 


1.08 


2.3 




—3.00 


GFGP 


3C 


67 


0.3012 


0.91 


0.99 


337.0 




—0.40 


vBMH 


3C 


79 


0.2559 


1.31 


1.01 


10.5 




—2.19 


SMP 


3C 


98 


0306 


4.97 


0.96 


9.0 




—2.75 


FB 


3C 


109 


3056 


1.64 


0.64 


320.0 




—0.71 


HI* 


4C 


14.11* 


0.2070 


0.89 


0.73 


29.7 




— 1.52 


HAPR 


3C 


123 


0.2177 


16.32 


0.83 


108.9 




—2.24 


HAPR 


3C 


132 


0.2140 


1.05 


1.00 


4.1 




—2.49 


HAPR 


3C 


153 


0.2769 


1.35 


0.87 


< 0.5 




—3.52 


HAPR 


3C 


171 


0.2384 


1.22 


0.80 


2.2 




—2.82 


GFGP 


3C 


173.1 


0.2920 


0.77 


1.01 


7.4 




—2.13 


GFGP 


3C 


184.1 


0.1182 


1.24 


0.73 


6.0 




—2.35 


RP 


DA 240* 


0350 


0.69 


0.95 


121.9 




—0.69 


MKOW 


3C 


192 


0598 


2.05 


0.71 


6.0 




—2.55 


H79 


4C 


14.27* 


0.3920 


0.39 


0.99 


< 0.5 




—3.03 


LP 


3C 


215 


0.4110 


0.41 


0.85 


20.0 




— 1.42 


PH 


3C 


219 


0.1744 


2.29 


1.04 


51.6 




— 1.71 


HAPR 


3C 


223 


0.1368 


1.29 


0.75 


9.0 




—2.20 


RP 


4C 


73.08* 


0.0581 


0.90 


0.84 


19.0 




— 1.69 


M79 


3C 


234 


0.1848 


1.54 


1.04 


90.0 




— 1.29 


RP 


3C 


244.1 


0.4280 


1.12 


0.89 


< 0.6 




—3.41 


GFGP 


3C 


249.1 


0.3110 


0.78 


0.94 


110.0 




-0.91 


PH 


3C 


268.3 


0.3710 


1.09 


0.96 


< 20.0 


< 


-1.86 


F90 


3C 


274.1 


0.4220 


0.76 


1.08 


< 6.0 


< 


-2.27 


JPR 


3C 


284 


0.2394 


0.69 


0.71 


3.2 




-2.40 


GFGP 


3C 


285 


0.0794 


0.76 


0.81 


6.8 




-2.07 


HAPR 


3C 


300 


0.2700 


1.10 


0.90 


9.0 




-2.18 


RP 


3C 


303 


0.1410 


0.94 


0.80 


150.0 




-0.78 


PH 


3C 


319 


0.1920 


0.65 


1.07 


< 1.0 


< 


-2.89 


BF 


3C 


321 


0.0960 


1.22 


0.83 


30.0 




-1.63 


JPR 


3C 


326 


0.0895 


0.48 


1.63 


3.5 




-2.20 


R90 


3C 


349 


0.2050 


1.14 


0.80 


24.0 




-1.73 


JPR 


3C 


351 


0.3710 


1.21 


0.84 


15.0 




-2.02 


L81 


3C 


381 


0.1605 


1.29 


0.94 


5.2 




-2.45 


F84 


3C 


3S2 


0.0578 


2.22 


0.74 


251.2 




-0.91 


HAPR 


3C 


388 


0.0908 


1.77 


0.91 


57.9 




-1.51 


HAPR 


3C 


390.3 


0.0569 


4.48 


0.62 


340.0 




-1.10 


vBJ 


3C 


401 


0.2010 


1.37 


1.14 


32.0 




-1.71 


BBDG 


3C 


436 


0.2145 


0.99 


0.99 


20.0 




-1.77 


RP 


3C 


438 


0.2900 


1.54 


1.20 


10.0 




-2.32 


L81 


3C 


452 


0.0811 


3.26 


0.95 


130.0 




-1.42 


RP 


3C 


457 


0.4270 


0.55 


0.93 


6.0 




-2.10 


LRL 




-4 -2 

log R 

Figure Al. The a posteriori probability density function of view- 
ing angle, 8, as a function of core-to-lobe ratio, R, under the 
assumption described in the text. The solid line indicates the 
maximum of the pdf, and the dashed lines show the lo~ limits. 



angle. Rather than adopt an ad hoc formulation to create 
a smooth transition at 8 « 9 C and prevent this unrealistic 
behaviour, we choose 6 C = 0. The behaviour o f OC R) under 
these assumptions is shown graphically in Fig. Al. 
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